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Abstract 

We present an evaluation of the quark mass renormalization factor for Nf = 2 + 1 QCD. 
The Schrodinger functional scheme is employed as the intermediate scheme to carry out non- 
perturbative running from the low energy region, where renormalization of bare mass is performed 
on the lattice, to deep in the high energy perturbative region, where the conversion to the renor- 
malization group invariant mass or the MS scheme is safely carried out. For numerical simulations 
we adopted the Iwasaki gauge action and non-perturbatively improved Wilson fermion action with 
the clover term. Seven renormalization scales are used to cover from low to high energy regions and 
three lattice spacings to take the continuum limit at each scale. The regularization independent 
step scaling function of the quark mass for the Nf = 2 + 1 QCD is obtained in the continuum 
limit. Renormalization factors for the pseudo scalar density and the axial vector current are also 
evaluated for the same action and the bare couplings as two recent large scale Nf = 2 + 1 simu- 
lations; previous work of the CP-PACS/JLQCD collaboration, which covered the up-down quark 
mass range heavier than m-jx ~ 500 MeV and that of PACS-CS collaboration for much lighter quark 
masses down to m-j^ = 155 MeV. The quark mass renormalization factor is used to renormalize 
bare PCAC masses in these simulations. 
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I. INTRODUCTION 



The strong coupling constant and the quark masses constitute the fundamental param- 
eters of the Standard Model. It is an important task of lattice QCD to determine these 
parameters using inputs at low energy scales such as hadron masses, meson decay con- 
stants and heavy quark potential quantities. The results can be compared with independent 
determinations from high energy experiments. 

In the course of evaluating these fundamental parameters we need the process of renor- 
malization in some scheme. The MS scheme is one of the most popular schemes, and hence 
one would like to evaluate the running coupling constant and quark masses through input 
of low energy quantities on the lattice and convert it to the MS scheme. A difficulty in 
this process is that the conversion is given only in a perturbative expansion, and should be 
performed at much higher energy scales than the QCD scale. At the same time the renor- 
malization scale /i should be kept much smaller than the lattice cut-off to reduce lattice 
artifacts. Furthermore we may need to suppress artifacts due to a finite lattice extension L. 
We therefore require 

T<AQCD</i<-. (I.l) 
L a 

The practical difficulty of satisfying these inequalities in numerical simulations is called the 
window problem. 

The Schrodinger functional (SF) scheme [l|-[9[ is designed to solve the window problem. 
It has an advantage that systematic errors can be unambiguously controlled. A unique 
renormalization scale is introduced through the box size L in the chiral limit and the scheme 
is mass independent. A wide range of renormalization scales can be covered by the step 
scaling function (SSF) technique. This matches our goal to obtain the coupling constant 
and quark masses in the MS scheme. The SF scheme has been applied for evaluation of the 
QCD coupling for = 0, iV^ = 2 ["i^J and = 2 + 1 @ cases. 

For the quark mass renormalization factor the SF scheme has been applied for Nf = jsj 
and Nf = 2 [sj QCD. At low energy scales of /i ~ 500 MeV, where physical input is given, 
we expect the strange quark contribution to be important in addition to those of the up 
and down quarks. Thus the aim of the present paper is to go one step further and evaluate 
the quark mass renormalization factor in Nf = 2 + 1 QCD. Our goal is to renormalize 
the bare light quark masses and meson decay constants evaluated in two recent large-scale 
Nf = 2 + 1 lattice QCD simulations employing non-perturbatively 0(a) improved Wilson 
quark action; the work of CP-PACS/JLQCD Collaboration with relatively heavy pion mass 
of ~ 500 MeV [l^ and that of PACS-CS collaboration with much lighter quark masses 



down to = 155 MeV 11, 12 



II. OUR STRATEGY 

Our goal is to evaluate the renormalization group invariant (RGI) quark mass M given 



by 
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where 'g{n) and m(yu) are renormalized coupling and quark mass in some scheme at scale 
/X. I3{g) and T{g) are the renormalization group (RG) function of coupling and mass in the 
same scheme. 

Once the RGI mass M, which is scheme independent, is evaluated, the conversion into 
the MS scheme can be carried out from the equation (lll.ip with the RG functions in the MS 
scheme. 

We therefore first derive the renormalization factor ZM{go), which converts the bare 
PCAC mass at a bare coupling go to the RGI mass. Since the step scaling function (SSF) 
cr(n) for the running coupling has already been obtained j9|, an evaluation of the RGI mass 
renormalization factor proceeds in the same way as given in Ref. js], [8|, which we briefly 
summarize here. 

1. We start by preparing the coupling step scaling function a{u), which gives the relation 
between the renormalized coupling constants when the renormalization scale is changed 
by a factor two, 

aiu) = \im(f{2L) (II.2) 



u=g^(L),m=0 



The renormalization scale is defined through a physical box size L where the renor- 
malized coupling u = g'^{L) is evaluated in the SF scheme. 

2. We define a reference scale Lmax through a renormalized coupling constant Uq = 
^^(Ljnax)- The value of uq is arbitrary as long as it is well in the low energy region 
to suppress lattice artifacts with a/Lmax ^ 1 and within the range covered by the 
SSF. We need some physical input measured separately in simulation at large volume 
to evaluate Lmax in physical units. In this paper we shall adopt the lattice spacing 
determined from hadron masses for this purpose. 

3. We then calculate the step scaling function (SSF) Ep for the pseudo scalar density, 
which is given by the ratio of two renormalization factors at different renormalization 
scales n = 1/L and 1/(2L) at the same bare coupling go 



a\ Zp{go,a/{2L)) 
T 



Lj Zp{go,a/L) 



(11.3) 

g-^ {L)=u,m=0 



where mass independent scheme in the massless limit is adopted. We used the same 
definition of the pseudo scalar density renormalization factor Zp as Ref. [sl, [sl . 

After taking the continuum limit with three lattice spacings 

MmEpfu,^) =ap{u) (II.4) 

and performing a polynomial fit for crp(u) we get the non-perturbative renormalization 
group flow of the quark mass. 

4. Since the SSF (yp{u) is given as a function of renormalized coupling u we need to 
prepare a sequence of couplings Ui in order to perform a non-perturbative running 
starting from uo = ^^(-Z^max)- Wj's differ by a factor two in the renormalization scale 
Lj = 2~*Lmax and is given by solvingthe SSF cr^Ui^i) = Ui. We find that we are well 
in the perturbative region for z > 5 [9(] . 
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5. Taking the product of the SSF ap{ui) at each scale we get a factor 

m(l/L„) 



= w 

m(l/Lmax) 



'-' = X{crp{u,) (II.5) 



i=l 



which represents the running of the renormahzed mass from Lmax to L„. We derive 
the renormahzation factor, which converts the quark mass in the SF scheme at the 
reference scale Lmax to the RGI mass 

M m(l/L„) M _i M 

-w„ =777^. (II.6) 



m(l/Lmax) m(l/Lmax) " m(l/L„)' 

The last factor M /rn{l/ Ln) is evaluated perturbatively for n > 5 with T{g) at two-loop 
order and (3{g) at three-loop order pjj] using 

6. The renormalized PCAC quark mass at scale /i = l/Lmax is given by 

m(l/Li^ax) = lim ^^^90,a/L) ^(bare)^^^^^_ 
a^O Zp(5(o,a/Lmax) 

The bare PCAC quark mass '^pcac(5'o) is obtained in some simulation at large volume 
carried out at a bare coupling qq. We need to calculate the renormahzation factor 
ZA{go,a/L) for the axial current jH, 23| and Zp{go,a/L^s^^) for the pseudo scalar 
density at the same bare coupling go in this simulation. The box size L for the axial 
current need not coincide with the box size Lmax for the pseudo scalar density since 
the former is scale independent in the continuum limit. We finally obtain 

Zm(c/0, a/^max) = =7TT? 7 ^7 77 7- (H-S) 

max J 

The factor Zji{gQ,a/L) is also used for renormahzation of meson decay constants. 

7. As a last step, we evaluate the mass renormahzation factor in MS scheme at a scale 

Z^igo, /i) = ^^ZMigo, a/L^^^) (II.9) 



The factor m^{n) / M is given by using fill. ID with the renormahzation group functions 
in MS scheme at four loops. The renormahzation scale n is introduced through the 



RGI scale Aj^jg evaluated in Ref . [9[ . 



III. NUMERICAL SETUPS 

We adopt the renormahzation group improved Iwasaki gauge action and the 0{a) im- 
proved Wilson fermion action with the clover term, whose parameters are identical to those 
adopted in the previous paper for the running coupling [oj : The boundary improvement co- 
efficients of the gauge action are set to the tree- level values cf = 1 and cf^ = 3/2 for Iwasaki 
action [isj . The improvement coefficient csw of the clover term is known non-perturbatively 



14i |. The boundary improvement coefficient q of the fermion action is set to the one loop 
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value [15|. Differences reside in the Dirichlet boundary condition for the spatial gauge link, 
which is set to 




(III.l) 



and the twisted periodic boundary condition of the fermion fields in the three spatial direc- 
tions 

^{x + Lk) = e'^i){x), ii{x + Lk) = e~'^^{x), (III.2) 

where 9 = 0.5 [sl, [1]. These conditions are known to show a good perturbative behavior [l^ 
for the quark mass renormalization factor. 

For numerical simulations, the HMC algorithm is adopted for up and down quarks and the 
RHMC algorithm for the strange quark. The masses of all three quarks are set to zero. We 
adopt the CPS-I-+ code [13] and add some modification for the SF formalism. Simulations 
were carried out mainly on the T2K-tsukuba computer at University of Tsukuba and the 
T2K-tokyo computer at University of Tokyo [18|] . 



IV. STEP SCALING FUNCTION 



A. Pseudo scalar density renormalization factor 

We adopt seven renormalized coupling values to cover weak (^^ = 1.001) to strong {g"^ = 
3.418) coupling regions [9)], which approximately satisfy gf_^i{L) = g]{2L). For each coupling 
we use three box sizes L/a = 4, 6, 8 to take the continuum limit. At the three lattice sizes the 
values of /3 and k were tuned to reproduce the same renormalized coupling keeping the PCAC 
mass to zero. On the same parameters we evaluate the pseudo scalar density renormalization 
factor Zp{go,a/L) and Zp{go,a/{2L)) with the Dirichlet boundary condition fllll.ip . 

We adopt the definition of the renormalization factor js], Isj given by 



Zp{go,a/L) = —^^^ (IV.l) 
np fp{L/2) 

where /p(xo) and /i are two-point functions of pseudo scalar density at bulk and at boundary 
given by 

fpM = -j^Y^ E {P%^)o-) , (IV.2) 

/i = -]^|^(^"^^")' (IV.3) 
P'^(a;)=V^(x)75^r"^(x). (IV.4) 

The boundary pseudo scalar operators O"" and O'"' are written in terms of the boundary 
quark fields C [4] as 

(^'^ = EC(^)75yC(^l, 0'^ = YX{u)lJ-^C\v). (IV.5) 
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The normalization factor 



Up = 1 + am + 31— cos — 



(IV.6) 



is chosen to give Zp{go = 0,a/L) = 1 at tree level on the lattice, where the quark mass am 
is written in terms of the PC AC mass as 



am 



1 — 31 — cos 



aO' 



\ 



2am 



PCAC 



'ia'^mpQj^^ + 6 sin^ 



aO 

T 



3 sin 



a9\ 
(IV.7) 



The normalization factor np is calculated by setting 9 = 0.5 and the PCAC mass at each 
(/?,«:). 

The value of the renormalization factors are listed in table [T] at each parameter together 
with the step scaling function (111.30 . In the table we notice that there is a tiny difference in 
the renormalized coupling g'^{L) between three boxes L/a = 4,6,8. We adopt the coupling 
at L/a = 8 to define the renormalization scale and the deviation is corrected perturbatively 



19|. Statistics of the runs are given in table HT] together with values of the PCAC mass 
defined by 



"^PCAC 



fA{xo) 



HVo + VS)/A(f)+CAVSVo/p(f) 



2ff 



imp. 



(IV.8) 
(IV.9) 



in terms of the improved axial vector current with non-perturbative improvement coefficient 
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21-0.195^2 



(IV.ll) 



B. Perturbative improvement and continuum limit 

We perform a perturbative improvement of the SSF before taking the continuum limit. 
For this purpose we need a perturbative evaluation of the lattice artifact in the SSF 

Sp{u, a/L) = ^pML)-ap{u) _ ^^^^^2) 
ap[u) 

For the SSF of the running coupling the artifact 6{u,a/L) was evaluated at one loop level 



13|, l21|. However, the magnitude of the one- loop correction was found to be large and hence 
the range of applicability is limited only for very weak coupling region for our lattice action, 
revealing the importance of two-loop coefficient [9i]. Since 6p is not known at all for our 
setup, instead of calculating 6p at one and two-loop level perturbatively we calculate SSF's 
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directly by Monte-Carlo simulations at very weak coupling /3 > 10. The same value of k is 
used as in Ref. [oj, which gives almost vanishing PCAC mass there. The results are listed in 
table Hm with the number of configurations and the value of PCAC mass given in table IIVI 
We define (5p(m, a/ L) by the deviation from the perturbative SSF's Up at two-loop order 
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5p{u,a/L) = ' ^ • (IV.13) 

ay[u) 

The deviation is fitted to a polynomial form for each a/L, 

1 + 5p{u, a/L) = 1 + di{a/L)u + d2{a/L)u^. (1V.14) 

We tried a quadratic fit using data at m < 1.524, which is plotted in Fig. [H Since the 
quadratic fit provides a reasonable description of data we opt to cancel the 0(a) contribution 
dividing out the SSF by the quadratic fit 

^ V' l) ~ l + d^{a/L)u + d2{a/L)u^- ^ ^ ""^ 

We notice that the deviation is consistent with zero within one standard deviation for 
5p(m, 1/8) even at m ^ 1. We consider a large deviation in the quadratic fit is an arti- 
fact due to relatively large statistical error for the data at m ^ 1 and we therefore do not 
apply the present improvement for a/L = 1/8. Even if we apply the improvement to the 
SSF 6p{u, 1/8), the quark mass renormalization factor changes only by (9(1%), which is 
within statistical errors. 

We give the values of the perturbatively improved SSF's in table |V] for three lattice 
spacings at each of the 7 renormalization scales. Scaling behavior of the improved SSF 
is plotted in Fig. [2l Almost no scaling violation is found. We performed three types of 
continuum extrapolation: a constant extrapolation with the finest two (filled symbols) or all 
three data points (open symbols), or a linear extrapolation with all three data points (open 
circles), which are consistent with each other. We employed the constant fit with the finest 
two data points to find our continuum value. 

The RG running of the continuum SSF is plotted in Fig. [31 A polynomial fit of the 
continuum SSF to third order yields 

ap{u) = 1 + pou + piu^ + P2U^, (IV.16) 
p^ = -0.002826, p2 = 0.000031 (IV.17) 

fixing the first coefficients to its perturbative value po = — 8 In 2/(4??)^. We notice that 
Pi is consistent with its perturbative value Pi^'^^ = —0.00282843 at two loops. The fitting 
function is also plotted (solid line) together with the two loops perturbative running (dashed 
line). 



C. Non-perturbative running mass 

Once m(l/Lmax)/^ is evaluated from the SSF given in flIV.161) according to the strat- 
egy in sectjni we are able to derive the non-perturbative running mass in the SF scheme 
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according to (1II.6|) . The running mass at scale /z = is given in units of the RGI mass 
M by 

m(l/L„) _ _-Jn{l/L 

max J /T-, r 1 r)\ 

M " " M ■ ^ ^ ^ 

The resuh is plotted in Fig. H] as a function of yu/AsF- The solid line is perturbative running 
given by using T{g) at two-loop order and P{g) at three- loop order in (lll.ll) with inputs 
of m{l/ Li2)M and g{l/Li2) at a very high energy scale of L12 = 2~^^Lmax- We omit the 
horizontal error bar, which may be introduced from that of Asf- 



V. Zp AT LOW ENERGY SCALE 

Some years ago CP-PACS and JLQCD Collaborations jointly performed an Nf = 2 + 1 
simulation with the 0{a) improved Wilson action and the Iwasaki gauge action 10|. Three 
values of f3, 1.83, 1.90 and 2.05 were adopted to take the continuum limit and the up-down 
quark mass covered a rather heavy region corresponding to m^/mp = 0.63 — 0.78. This 
project has been taken over by the PACS-CS Collaboration aiming at simulations at the 
physical light quark masses [III, At present result at a single lattice spacing /3 = 1.90 is 



available with very light quark masses down to niTr/mp ^ 0.2. 

One of the purpose of this paper is to provide the non-perturbative renormalization factor 
to renormalize the bare PCAC quark masses measured in these large scale simulations. We 
start by introducing a reference scale Lmax- For this purpose, the renormalized coupling 
is evaluated at the same bare coupling /3 as in the CP-PACS/ JLQCD simulations and in 
the chiral limit j9|. The reference scale Lmax is given by the box size we adopt in this 
evaluation. The renormalized coupling ^^(Lmax) should not exceed our maximal value 5.13 
for the coupling SSF significantly. We use the box size of L/a = 4 for /3 = 1.83 and 1.90 to 
define L^^x and L/a = 6 for /3 = 2.05. A physical scale is introduced into the present work 
so that the reference scale is translated into MeV units. We employ three types of hadron 



masses m^, irip, uik (^^0) uM "^w, mx, [III, ll2[ as inputs and use the lattice spacing 
a as an intermediate scale, which is given in Table IVII 

Then we evaluate the pseudo scalar density renormalization factor at the same /3 = G/g^ 
and K as that in Ref. fsj to give the denominator of the PCAC mass renormalization factor 
(]II.7p . The values of the renormalized coupling constant j9| and the renormalization factor 
Zp{gQ,a/L^ax) are listed in Table rvTTl 



VI. AXIAL CURRENT RENORMALIZATION FACTOR 



The remaining ingredient of the renormalization factor ( 111.71) at low energy is that of the 
axial vector current. We calculate the renormalization factor according to the procedure 



in Ref. [22 
condition 



23 



23( 1 through the axial Ward-Takahashi identity. We adopt the renormalization 



ZAigo, m) = Za[ go, m, -T 
Za (go, m, xq) 



fi 



"-^ Jaa {^0, 1^) - 2mpcAc/i'A {^o, 1^ 



(VI.l) 
(VI.2) 
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which is apphcable to small non-vanishing PCAC mass. In this section we suppress the 
a/-Z^max dependence while we explicitly write the m dependence instead. Here /^^ and /p^ 
are four-point functions of boundary operators and the improved axial current (IIV.IOI) or 
the pseudo scalar density operator (IIV.4p defined by 

/iA(^o, 2/0) = E ^"""^^'^ (o" {^rY (x) {j^r)\y)o^ , (VI.3) 

xy 

fU^o, yo) = E e'^'^e^'' (o'^P'^ix) {Ar-y iy)0^ . (VI.4) 

xy 

The PCAC mass here is given by an average over three time slices 

1 tt' -2 (Vo + V;) IaM + CAV'„Vofp(xo) 

The normalization factor is chosen to give Za^Qo = 0,m, Xq) = 1 at tree level on the 
lattice. 

As a byproduct we are able to derive the renormalization factor for the vector current 

V,{x)=ii{x)^,Ky{x) (VI.6) 

with negligible computational cost. Here again we adopt the same definition for the renor- 
malization factor as in [22[ 23 1 given by 



Zy{g,) = Zy (^go, I) , (VI.7) 

Zv{go,Xo) = —Y^, (VI.8) 

Mxo) = ^ E '^'^^ {o'Xi^)<^') ■ (vi-9) 

X 

The normalization factor is given by ny = np. 

There are several alternative choices in the calculation of Z^. box size L/a, twisting 

angle 6 and two definitions of Za_ in Ref. H, 23 1 with or without the disconnected diagrams. 



Our requirement is that the 0(0^) lattice artifact is suppressed as much as possible within 
a reasonable computational cost. As a measure of the lattice artifact, we use the time 
dependence of Za/V: which should disappear in the continuum limit. 

We have varied the spatial box size from L/a = 6 to 12 until we find a plateau in Z^/v 
as a function of Xq. Our box sizes are 6^ x 18, 8^ x 18, 10'^ x 24 and 12^ x 30. The temporal 
lattice size was taken to be a multiple of six so that it is even and the two axial currents 
can be placed at an equal distance T/ ( 3a) from the boundary and from each other. We set 
T/a to the nearest value to 9/4 x L/a [23|. The spatial twist angle is set to 6' = and 0.5 
to examine its dependence. The bare quark mass is tuned so that |ampcAc| < 0.01 since 
the dependence is known to be negligible [23[ in the definition (1VI.2P for such a small mass. 



The Dirichlet boundary condition f llll.ip is the same as that for Zp. 

The time dependence of ZA{go,m,xo) is plotted for (3 = 1.83 in Figs. [5] (6' = 0.5) and [6] 
= 0). We found no plateau for both definitions with or without the disconnected diagram 
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for almost all box sizes. The variation of 9 does not change the conclusion. An exception is 
that with the disconnected diagram on a 12^ x 30 lattice indicated by red diamonds. However 
we cannot be very confident of a plateau for two reasons. One is that the statistical error is 
too large and the other is a large discrepancy between the two definitions with or without 
the disconnected diagrams, which should be a pure lattice artifact at the massless point. 

The time dependence is more tamed at /3 = 1.90 as is seen in Figs. [7] and [SI We found 
a plateau in with the disconnected diagram for L/a > 10 at ^ = 0.5 or L/a > 8 with 
^ = 0. The value of Za is consistent with each other at = 2T/3 within one standard 
deviation once a plateau is located. We have no trouble to find a plateau for /3 > 2.05 as 
is shown in Figs. [9l [TOj Furthermore results from the two definitions with or without the 
disconnected diagrams are consistent with each other. 



We conclude that we should adopt the definition with the disconnected diagrams [23 
to locate a plateau at smaller /3. It is better to set the twist angle 9 to zero although 
computational cost is slightly heavier and statistical fluctuation becomes worse. The box 
size should be taken as large as possible since lattice artifacts should manifest as a polynomial 
of a/L. 

In this paper we adopted the following box size at each /3 to define Za'- 12^ x 30, 9 = 0.5 
at /3 = 1.83, 10^ X 24, = at /3 = 1.90 and 12^ x30,9 = 0.5 for (3 = 2.05. The disconnected 
contribution is included. The physical box size varies from 0.84 fm to 1.3 fm. We derive 
our Za by fitting the plateau around xq = 2T/3 by a constant. The fitting range is plotted 
in Figs. [5] - [TOl The error is estimated by the Jackknife method. The results are listed in 
tables IVIIIl HXl and |X] together with those at various box sizes. 

In Fig. [TT]we plot the values for Za according to our definition (filled circles), together 
with those of other definitions (open symbols). Scattering of points starting around /3 = 1.95 
indicates that lattice artifacts are increasingly large in our data for large lattice spacings. 

The time dependence of Zv{go,xo) is plotted for (3 = 1.83 in Fig. [T2]for 9 = 0.5 (left) and 
9 = (right) with horizontal axis normalized to unity. For L/a> 10 we find a long plateau 
whose values does not depend on L and 9. The behavior is almost the same for /3 = 1.90 
(Fig. [T3D and /3 = 2.05 (Fig. [H]). Zyigo, xo) seems to be almost fiat at /3 > 3.0 for L/a = 8 
as is seen in Fig. [T3 The renormalization factor Zy given by fitting the plateau is listed in 
tables Emi EX] and H 

In Fig. [in] we plot our final results for ZA^go) and Zv{go) as a function of (3 together 
with perturbative behavior (solid line) and results from the tadpole improved perturbation 
theory (star symbols). 



VII. RGI MASS RENORMALIZATION FACTOR 

We derive the renormalization factor Zm for the RGI mass according to the procedure 
given in Sec. [TT] This factor is intended to renormalize the bare PCAC masses obtained in the 
two large scale simulations carried out at three values of /3. The three hadron masses, m.^, 
f^K, t^n were used in Ref. [U, 13] to determine the hght quark masses and the lattice spacing. 



Two choices m^r, rup ttik or m^r, rrip m(i, were adopted in Ref. [10[. The results for Zm are 
listed in Table [Xll Also listed are the mass renormalization factors Z^^(/3, /i = 2 GeV) in the 
MS scheme at a renormalization scale fi = 2 GeV. We emphasize that the renormalization 
factor here is defined in terms of the renormalization group functions for three fiavors. 
Results from tadpole improved perturbation theory is also listed for comparison. 
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The error in the renormahzation factor includes all the statistical and systematic ones 
except for that from the choice of the reference scale Lmax- Lattice artifacts may be present 
at order C((a/Lmax)^)- In order to estimate this effect we increase Lmax/o at each (3. There 
is a problem that values of the renormalized coupling Wmax and the renormahzation factor 
'^p(-f'max) tend to exceed the region covered by our step scaling function as is shown in 
Table IXIII Since we attempt just a rough estimate, we adopted the step scaling function 
by extrapolating its polynomial form. The result is listed in Table IXIIII and by comparing 
with values in Table IXll we find a few percent effect at /3 = 2.05, while the magnitude may 
increase to a 10 % level at lower values of /3. 

As the last step we apply our renormahzation factor to the bare PCAC masses in Refs. [13- 



12[ | to obtain values for the renormalized quark masses. The numerical results are given in 



Table IXIVI and are plotted in Fig. [T7] both for the averaged up and down quark mass (left 
panel) and for the strange quark mass (right panel). For the old CP-PACS/JLQCD work 
of Ref. [lOf results with K and (j) meson input are plotted (filled squares and triangles) 
together with perturbatively renormalized masses using tadpole improved renormahzation 
factor (open squares and circles). The triangle represents the result for the more recent work 
of PACS-CS with the pion mass reaching down to m,r = 155 MeV. 

It is disappointing that the old CP-PACS/JLQCD results do not exhibit a better scaling 
behavior by going from perturbative to non-perturbative renormahzation factor. However, 
we should note a significant change in the average up and down quark mass with the recent 
PACS-CS work (filled triangle). This represents a systematic error due to chiral extrapola- 
tion of the old CP-PACS/JLQCD work whose pion mass reached only ~ 500 MeV. We 
should also note that the renormahzation factor Za involves a large uncertainty at /3 = 1.83 
which is not reflected in the error bar of Fig. [T71 We feel that results at /3 = 2.05 using 
simulation with physical pion mass are needed to find the convincing values for light quark 
masses with our approach. 



VIII. CONCLUSION 

We have presented a calculation of the quark mass renormahzation factor for the Nf = 
2 + 1 QCD in the mass independent Schrodinger functional scheme in the chiral limit. We 
adopt seven renormahzation scales to cover from low to high energy region. Three lattice 
spacings are used to take the continuum limit. We calculate the pseudo scalar density 
renormahzation factor and obtain the SSF for the running mass. Applying a perturbative 
improvement, we find that the step scaling function shows a good scaling behavior and the 
continuum limit may be taken safely with a constant extrapolation of the results for the 
finest two lattice spacings. The non-perturbative SSF turned out to be almost consistent 
with the perturbative two loops result. 

With the non-perturbative renormahzation group flow we are able to estimate the quark 
mass renormahzation factor which renormalizes the bare PCAC mass and give the RGI 
mass. For this purpose we derive the renormahzation factor of the pseudo scalar density 
at low energy and that of axial vector current. The results are used to renormalized the 



PCAC masses in the two recent large scale simulations |10l4l2|. We also evaluate the non- 
perturbative factor Z^^{f3,fi = 2 GeV) in MS scheme at a scale /i = 2 GeV. 

All statistical and systematic errors are included in the renormahzation factor Zm except 
for that due to the choice of the reference scale L^ax- A rough estimate of C((a/Lmax)^) 
effect suggests possible 10 % systematic error at /3 = 1.90. However, a firmer conclusion 
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requires the step scaling function at the couphngs stronger than those explored in the present 
work. We leave it as future work to examine this issue. 

Applying our non-perturbative renormalization factor to the present PACS-CS simulation 
with the pion mass as light as m-Tr = 155 MeV yields m^^ ifi = 2 GeV) = 2.78(27) MeV 
for the average up and down quark, and m^^(/i = 2 GeV) = 86.7(2.3) MeV for the strange 
quark. Simulations at weaker couplings under way will tell if these values stay toward the 
continuum limit. 
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TABLE I: The pseudo scalar density renormalization factor Zp at scale L and 2L. The value of 
f5 and k has been tuned to reproduce the same physical box size L and near zero PCAC mass in 
Ref. The step scaling function Sp is also listed. 
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TABLE II: Number of configurations for each run. Also listed is the PCAC mass at each parameter. 



/3 


K 


L/a 


u = tiL) 


Zp{go,L/a) 


2L/a 


Zp{go,2L/a) 


Sp(u, a/L) 


10 


0.1270893 


4 


0.58565(34) 


0.95007(21) 


8 


0.92391(27) 


0.97247(36) 


20 


0.1260654 


4 


0.29543(17) 


0.975596(93) 


8 


0.96211(11) 


0.98618(15) 


40 


0.1255571 


4 


0.14876(23) 


0.987683(80) 


8 


0.981223(64) 


0.99346(10) 


60 


0.1253871 


4 


0.099336(40) 


0.991832(46) 


8 


0.9873514(42) 


0.995482(62) 


10 


0.1272305 


6 


0.59707(44) 


0.93262(13) 


12 


0.91141(31) 


0.97726(36) 


20 


0.1261216 


6 


0.29775(33) 


0.96681(11) 


12 


0.95557(13) 


0.98837(18) 


40 


0.1255700 


6 


0.149219(45) 


0.983422(56) 


12 


0.977882(57) 


0.994366(81) 


60 


0.1253863 


6 


0.099664(47) 


0.988962(38) 


12 


0.985302(42) 


0.996299(57) 


10 


0.1272310 


8 


0.6051(12) 


0.92377(19) 


16 


0.90192(77) 


0.97635(85) 


20 


0.1261176 


8 


0.29971(49) 


0.962064(82) 


16 


0.95088(29) 


0.98837(31) 


40 


0.1255626 


8 


0.14948(27) 


0.981092(50) 


16 


0.97552(10) 


0.99432(12) 



TABLE III: The pseudo scalar density renormalization factor Zp at scale L and 2L at very weak 
coupling region /? > 10. The step scaling function Sp is also listed. 
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TABLE IV: Number of configurations for each run at very weak coupling. Also listed is the PCAC 
mass at each parameter. 
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TABLE V: The perturbatively improved SSF T,p {u,a/L) at "two loop" level. The improvement 
is not applied for a/L = 1/8. The SSF ap{u) in the continuum is also listed, which is given by a 
constant fit of two data at finest lattice spacings 1/6, 1/8. 
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TABLE VI: Lattice spacing a in unit of fm from large scale simulations [10| and [Hi, [1 
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TABLE VII: The renormalized coupling and the renormahzation factor Zp at ^ = 1.83, 1.90, 2.05 
for low 
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TABLE VIII: The (axial) vector current renormahzation factor at /? = 1.83 for various box sizes. 
Two choices oi Q = 0.5 and ^ = are adopted for comparison. We also listed Za without discon- 
nected diagrams as Z^""^. Data with star symbol is adopted for our Za- 
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TABLE IX: The (axial) vector current renormahzation factor at /3 = 1.90 for various box sizes. Two 
choices of f? = 0.5 and ^? = are adopted for comparison. We also listed Za without disconnected 
diagrams as Z^""^. Data with star symbol is adopted for our Za- 
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0.93183(39) 


0.93154(31) 



TABLE X: The (axial) vector current renormahzation factor at /? = 2.05 and larger for 9 = 0.5. 
We also listed Za without disconnected diagrams as Data with star symbol is adopted for 

our Za- 
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p 






Fi(taclj Z„ (K) 


Zm W) 




1.83 [10] 
1.90 [10] 
2.05 [10] 


1.33(15) 
1.693(46) 
1.862(41) 


1.04(12) 
1.315(35) 
1.417(29) 


1.07161 
1.09973 
1.14487 


1.04(12) 
1.315(35) 
1.416(29) 


1.07019 
1.09955 
1.1446 


1.90 [12] 


1.693(46) 


1.347(36) 


1.11322 







TABLE XI: Zm for the RGI mass and Z^^{2 GeV) in the MS scheme, which are designed to 
renormahze the bare PCAC mass in two recent simulations of Ref. 10| and (K) or 

(i;^) means which meson mass is used for physical scale input. Perturbative results with tadpole 
improvement are also hsted for Z^^{2 GeV). A slight difference in Z^^ at /3 = 1.90 between 
Ref. [10| and [Ij] is due to that of the lattice spacing, from which physical unit is introduced. 





K 


Lma.x/ 0- 




W-PCAC 


Zp 


"IpcAC 


1.83 


0.138192 


6 


9.37(27) 


0.00231(59) 


0.44601(86) 


0.00439(35) 


1.90 


0.137289 


6 


6.60(12) 


0.00035(37) 


0.49605(66) 


0.00313(22) 


2.05 


0.13611599 


8 


6.01(21) 


0.00026(15) 


0.51155(57) 


0.000733(74) 



TABLE XII: The renormalized coupling and the renormalization factor Zp at larger box size of 
Lmax/a at three /3's. Two PCAC masses are evaluated with different Dirichlet boundary condition. 
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Zm 


(K) 


1.83 [lOj 
1.90 [10] 
2.05 [10] 


1.46(18) 
1.837(57) 
1.890(47) 


1.15(14) 
1.418(41) 
1.437(31) 


1.90 [12] 


1.837(57) 


1.452(42) 



TABLE XIII: Zm and Z^^{2 GeV) evaluated with larger box size L, 



/3 


M^^' (K) 


m^i {K) 


Mf^^ (K) 


mf^ {K) 




{<t>) 


Mf^^ {cp) 


mf^ ((/)) 


1.83 [10] 
1.90 [10] 
2.05 [10] 


3.30(38) 
4.47(17) 
5.29(24) 


2.59(30) 
3.47(13) 
4.02(18) 


86(10) 
115.3(4.4) 
136.1(6.5) 


67.7(7.9) 
89.5(3.4) 
103.6(4.9) 


3.31(39) 
4.46(16) 
5.29(24) 


2.58(30) 
3.46(12) 
4.02(18) 


99(11) 
124.8(5.4) 
143.0(8.1) 


77.7(9.0) 
96.9(4.2) 
108.8(6.1) 


1.90 [12] 


3.49(34) 


2.78(27) 


109.0(3.0) 


86.7(2.3) 











TABLE XIV: Non-perturbatively renormalized mass for the averaged up and down quark and for 
the strange quark in the two recent simulations of Ref. 10l4l2l] . M^^^ is the RGI mass and is 



that in the MS scheme at a scale /x = 2 GeV. The unit is in MeV. {K) or {(f)) means which meson 
mass is used for physical scale input. 
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FIG. 1: The polynomial fit of the discrepancy Sp {u,a/L) /ap (u) at high /3 ^ 4. The fit is given 
for each lattice spacings a/L = 1/4 (left), 1/6 (middle) and 1/8 (right). Solid line is a quadratic 
fit. 
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FIG. 2: The SSF on the lattice with its continuum extrapolation at each renormalization scale. 
Three types of continuum extrapolation are tried: a constant extrapolation with the finest two 
(filled symbols) or all three data points (open symbols), or a linear extrapolation with all three 
data points (open circles). These are consistent with each other. 
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FIG. 3: The RG flow of the SSF. Dotted line is two loops perturbative running. Solid line is 
polynomial fit of the SSF. 
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FIG. 4: The non-perturbative running mass in the SF scheme. Solid line is a perturbative runnin 
with two and three loops RG function for the mass and the coupling. 
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FIG. 5: xq dependence of ZA{go,m,xo) at /3 = 1.83 and 9 = 0.5 for box sizes 6^ X 18 (upper left), 
8^ X 18 (upper middle), 10^ x 24 (upper right) and 12^ x 30 (lower). Red diamonds and black 
circles are definition with or without the disconnected diagrams. Solid and dashed lines represent 
an expected plateau, which is given by fitting data around xq = 2T/3 by a constant. Error is 
estimated with the Jackknife method. 
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FIG. 6: Xq dependence of ZA{go,m, xq) at /3 = 1.83 and ^ = for box sizes 6^ X 18 (left), 8^ X 18 
(middle) and 10^ x 24 (right). 
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FIG. 7: xq dependence of ZAigo^rn^ xq) at /3 = 1.90 and 9 = 0.5 for box sizes 8^x18 (left), 10^x24 
(middle) and 12'^ x 30 (right). Red diamonds and black circles are definition with or without the 
disconnected diagrams. Solid and dashed lines represent an expected plateau, which is given by 
fitting data around xq = 2T/3 by a constant. Error is estimated with the Jackknife method. 
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FIG. 8: Xq dependence of ZA{go,'rn,XQ) at /? = 1.90 and ^ = for box sizes 8'^ x 18 (left) and 
10^ X 24 (right). 
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FIG. 9: Xq dependence of ZA{go,rn,XQ) at 13 = 2.05 and 9 = 0.5 for box sizes 8'^ x 18 (left) and 
12^ X 30 (right). 
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FIG. 10: xq dependence of ^^(^0,™-, xq) at /3 = 3.0 (left), 4.0 (middle) and 5.0 (right). The box 
size is 8'^ X 18 and twist angle is set to 6* = 0.5. 
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FIG. 11: /3 dependence of ZA{go) with disconnected diagrams (filled circles) and that without them 
(open circles). 
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FIG. 12: xo dependence of Zy(g(o,xo) at /3 = 1.83, 6 = 0.5 (left) and 6 = (right) for box 
sizes 6^ X 18 (green triangles), 8^ x 18 (blue diamonds), 10^ x 24 (red squares) and 12^^ x 30 (black 
circles). Solid and dashed lines represent an expected plateau, which is given by fitting data around 
xq = T/2 by a constant. 
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FIG. 13: Xq dependence of ZY{go^xo) at /3 = 1.90, = 0.5 (left) and 6 = (right) for box sizes 
8^ X 18 (blue diamonds), 10^ x 24 (red squares) and 12^ x 30 (black circles). 
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FIG. 14: Xq dependence of Zy{gQ,m,XQ) at /3 = 2.05 and 9 = 0.5 for box sizes 8'^ x 18 (black 
squares) and 12^ x 30 (red circles). 
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FIG. 15: xq dependence of Zy(^o?^o) at ^ = 3.0 (left), /3 
box size 8^ x 18 and a value of 9 = 0.5 is adopted. 
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FIG. 16: P dependence of ZA{go) (left) and Zy{gQ) (right). Red filled circles are our non- 
perturbative renormalization factor. Solid line is a perturbative result at one loop. Stars are 
from tadpole improved perturbation theory. 
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FIG. 17: Scaling behavior of m^^ (left) and m^^ (right) 
PACS/JLQCD result 0, filled triangle is that of PACS-CS 
renormalized masses in Ref. 0. 
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